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Abstract. This paper considers the slow motion of the shock layer exhibited by 
the solution to the initial-boundary value problem for a scalar hyperbolic system 
with relaxation. Such behavior, known as metastable dynamics, is related to the 
presence of a first small eigenvalue for the linearized operator around an equilib- 
rium state; as a consequence, the time-dependent solution approaches its steady 
state in an asymptotically exponentially long time interval. In this contest, both 
rigorous and asymptotic approaches are used to analyze such slow motion for the 
Jin-Xin system. To describe this dynamics we derive an ODE for the the position 
of the internal transition layer, proving how it drifts towards the equilibrium lo- 
cation with a speed rate that is exponentially slow. These analytical results are 
also validated by numerical computations. 

Key words. Metastability, slow motion, internal layers, relaxation systems. 



The slow motion of internal shock layers has been recently widely studied. Such 
phenomenon, known as metastability, is usually related to the presence of a first 
small eigenvalue for the linearized operator around a given equilibrium state. From 
a general point of view, a metastable behavior appears when solutions exhibit a first 
time scale in which they are close to some non-stationary state for an exponentially 
long time before converging to their asymptotic limit. As a consequence, two dif- 
ferent time scales emerge: a first transient phase where a pattern of internal shock 
layers is formed in a 0{1) time scale, and a subsequent exponentially slow motion 
where the layers drift toward their asymptotic limits. 

A large class of evolutive PDEs, concerning many different areas, exhibits this 
behavior. Among others, we include viscous shock problems (see [13], [21], [25]). 
phase transition problems described by the AUen-Cahn equation, with the funda- 
mental contributions [5], [9], and Cahn-Hilliard equation, studied in [1] and [23] . 

We want to study the slow motion of the shock layer for the scalar hyperbolic 
system with relaxation 
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where the space variable x belongs to a one-dimensional interval / = {—£,£), i > 0. 
System (1.1) is a particular case of a class of more general hyperbolic relaxation 
systems of the form 



a 



u 



g{u,v) 



e ^q{u,v) 



usually utilized to model a variety of non equilibrium processes in continuum me- 
chanics: for example, non-thermal equilibrium gas dynamics ([15], [I9]), traffic dy- 
namics (|2], [IB], [IE]), and multiphase flows ([3], jl], Here e is a positive 
parameter, usually small, determining relaxation time. 

In the case of system (1.1), the parameter e can be seen as a viscosity coefficient] 
we are interested in studying the behavior of the solution to (1.1) in the limit of 
small e, and we want to identify the role of this parameter in the appearance and/or 
disappearance of phenomena of metastability. 

The main example we have in mind is the initial-boundary value problem for 
the quasilinear Jin-Xin system in the bounded interval / = (—£,£), with Dirichlet 
boundary conditions, that is 

dtU + dxV = 

1 



(1.2) 

for some £, 
(1.3) 



dtv + a^dxU = ^{f{u) 
u{±i, t) = u± 

u{x, 0) = uo{x), v{x, 0) = Vo{x) = /(mo(x)) 
a > 0, M-t G M and flux function / that satisfies 
f"is)>Co>0, /'(n+) < < /'(n_), /(«+): 



xe I, t>0 

t > 
X e I 

/('"-) 



We stress that, once the boundary conditions for the function u are chosen, the 
boundary conditions for the function v are univocally determined. This model was 
firstly introduced in [TU] as a numerical scheme approximating solutions of the hy- 
perbohc conservation law dtu + dxf{u) = 0. System (1.2) is strictly hyperbohc, with 
the spectrum a{df, dgY composed by two distinct real eigenvalues ±a. 

In the relaxation limit (e — )■ 0^), system (1.2) can be approximated to leading 
order by 

dtu + dxV = 



(1.4) 

that is 
(1.5) 

together with v 
(1.6) 



fin) 



dtU + dxf{u) = 
and complemented with boundary conditions 
u{—i,t)=u^ and u{£,t) = u+ 



SLOW MOTION OF INTERNAL SHOCK LAYERS FOR THE JIN-XIN SYSTEM 



3 



From the standard theory of entropy solutions to first-order quasi-hnear equations 



of hyperbohc type, it is known that (1.5) admits a class of solutions, hence possible 



discontinuous, with speed of propagation s given by the Rankine-Hugoniot condition 

^ _ 1/(^)1 
1^1 



Assumptions (1.3) guarantee that the jump from the value m_ to the value m+ is 
admissible if and only if «_ > and its speed of propagation s is equal to zero. 



In this case, equation (1.5) admits a large class of stationary solutions satisfying the 



boundary conditions, given by all that piecewise constant functions in the form 

X e (xo, i) 



uix] 



where xo is a certain point in the interval. Hence, given ^ G {—£,£), we can con- 
struct a one-parameter family {t^hyp(-;0} of steady states, parametrized by ^ that 
represents the location of the jump, and given by 

where xi denotes the characteristic function of the interval /. We remark that, 
once f/hyp(-; is chosen, the class of stationary solutions (f/hyp, ^hyp) for the original 
system (1.4) is given by the relation \4yp = /(f^hyp), so that 

For the initial-boundary value problem (1.5)-([L6), it is possible to prove that every 
entropy solution converges in finite time to an element of the family {t^hyp(";0}- 
This stabilization propriety has been proved for the first time in [l7], using the 
theory of generalized characteristic, firstly introduced in [6j. In this framework, 
assumptions (1.3) on the flux function / are crucial (see pUl Theorem 6.1]). Hence, 
every entropy solution to the initial-boundary value problem 

dtu + d^f{u) = 0, V = f{u) 
u{±i, t) = u± 

converges in finite time to an element of the family {Uhyp{-;^), Vhyp(-; 0}- 

For e > 0, the situation is very different. If we differentiate with respect to x the 



second equation of (1.2), we obtain 

(1.7) ut = e{a'^da;xU - duu) - d^f{u) 
Thus, stationary solutions to ( |1.2 ) solve 

(1.8) o?ed^^u = d^f{u) 

together with d^v = 0. The presence of the Laplace operator has the effect that 
only a single stationary state is admitted (see [12]). As an example, we consider the 
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case of Burgers flux, i.e. a = 1, f{u) = i^u . We can explicitly write the stationary 



solution for the problem (1.8)-(1.6) as 



(1.9) t/L.W = -fctanh(^) 

where k = k{e,i,u±) is implicitly defined by imposing the boundary conditions 
Moreover, V^^ri^) is defined by V^uri^) — /(^)- 

In the limit e — )■ 0"*", the single steady state (Ubury^bur) converges pointwise to 



{Uhjp{-', 0), Vhyp(-; 0)), while, for a class of general f{u) that verify hypotheses (1.3), 
the stationary solution {U'^,V'^) converges pointwise to (t/hyp(-; 0; ^yp('j 0); 
some ^ G /. 

Finally, the single steady state (f/^, V"^) is asymptotically stable (for more details 
see the spectral analysis performed in Section [s]), i.e. starting from an initial datum 
close to the equilibrium configuration, the time dependent solution approaches the 
steady state for t — )■ +oo. 

Next question is what happens to the dynamics generated by an initial datum 
localized far from the equilibrium solution (£/^, V^). Numerical computations show 
that, starting, for example, with a decreasing initial datum Uo{x) (see Figjl]), because 
of the viscosity, a shock layer is formed in a 0{1) time scale. More precisely, the 
solution generated by such initial datum still presents a smooth transition from m_ to 
u+, but the shock is located far away from zero, so that the solution is approximately 
given by a translation of the (unique) stationary solution of the problem. Once the 
shock layer is formed, it moves towards the equilibrium solution, and this motion 
is exponentially slow. Thus we have a first transient phase where the shock layer is 
formed, and an exponentially long time interval where the shock layer approaches 
the equilibrium solution. 

Concerning the function v, starting with the initial datum Vo{x) = f{uo{x)), 
we can observe that the position of the shock of u corresponds to the location 
the minimum value of the function v; so we have a first transient phase in which 
the profile of v stabilizes, and an exponentially slow phase where the value of the 
minimum of such profile drifts towards the value ^ that represents the location of 
the equilibrium solution for u. 

The aim of this paper is to study the dynamics generated by an initial datum 
localized far from the equilibrium solution and to determine a detailed description 
of the low-viscosity behavior of the solutions. 



To the best of our knowledge, the problem of the slow motion for the hyperbolic- 
parabolic Jin-Xin system (1.1) has been never investigated before. However, system 
(1.2) can be reduced by differentiation to (1.7) together with the equation dtU+dxV = 
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Figure 1. Profiles of {u,v), solutions to ([O]), with f{u) = a 
=1 e = 0.04 and u± = =pl. The initial datum is given by the couple 
(no(x), /(no(x))), with uo{x) a decreasing function connecting u+ and u-. 
Profiles at times t = 0, 0.2, 10, 10^ 0.5 x 10^. 



0, and, as stressed, the study of stationary solutions to (1.7) is the same of that of 
the scalar conservation law 



:i.io) 



dtu + dxf{u) = ed^^u 



together with the additional condition d^v = 0. 

A pioneering article that analyzes the dynamics of (1.10) for initial data close 
to the equilibrium solution has been published by G. Kreiss and H.-O. Kreiss \12\ . 
Here the flux function / is given by f{u) = so that equation (1.10) becomes 

the so-called viscous Burgers equation. To study the dynamics generated by such 
initial configurations, the authors consider the linearized equation close to the unique 
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stationary state U^^^ defined in (1.9) , that is 



dtu = CeU := ed^xU + dx{a{x)u), a{x) := - f (U^uri^)) 

In [12] it is shown that, if f{u-) = the eigenvalues of Cs are all real and 

negative. Moreover 

= 0{e-^/') and < -c/e < V A; > 2 

This precise distribution shows that the large time behavior of solutions is described 
by terms of order e'^i*, so that the convergence to the asymptotically stable state 
U^^^ is exponentially slow, when e is small. 



Concerning the phenomenon of metastability for equation (1.10), such problem 
has been examined, among others, in [21] and in [H] . Here the different approaches 
are based either on projection method or on WKB expansion, but the common aim 
is to derive an equation for the position of the shock layer ^, considered as a function 
of time, that describes its slow motion towards the equilibrium location. In both 
papers, the analysis is carried on at a formal level and numerically validated. 

A rigorous analysis has been performed firstly in [7] (and generalized to the case 
of non-convex fiux in [S]). There, to study the slow motion of the internal layer, a 
one-parameter family of functions that approximate the stationary solution is chosen 
family of traveling waves with small velocity. 



The phenomenon of metastability for the equation (1.10) has been analyzed by 
C. Mascia-M. Strani in [20]. To study the problem of the slow motion, the au- 
thors introduce a one-parameter family of functions ^)}^g/, approximating a 
stationary solution U^{x). 

Hence, considering the linearized equation around f/^, it is shown that the eigen- 
values of the linearized operator verify 

(1.11) - Ce"'' < < 0, A|. < -C/e VA; > 2 

Moreover, the position of the shock layer ^(t) satisfies \^{t) — ^| < l^ok"'^ *, where 
~ e~^f^ . This estimate shows that the shock layer drifts toward the equilibrium 

solution with a speed rate proportional to the first eigenvalue A^, so that this motion 

is exponentially slow. 

Motivated by the analogies among the study of our problem and some results for 



the scalar conservation law (1.10), in this article we follow the approach presented 
in [2ni. 

- We build-up a one parameter family of approximate steady states 

{w^(x;OW = {f/^(-;0,v^^(-;e)W 

such that {U^{-] ^), V^^(-, ^)) := {U^ , V'^) for some ^, and with the additional propriety 
that (?7'^(-; ^), V^^(-, ^)) — )■ (t^hyp(-; 0^ ^hyp(-; 0) as e — )■ in an appropriate sense. 
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Moreover we require the error 

'i[W^]\ 

to be small in e in a sense to be specified. 

- We describe the dynamics of the system in a neighborhood of the family 

{f/^(-;0,^^(-;0}- 

Once a set of reference states 0^ ^^('S 0) is chosen, we determine spectral 

proprieties of the linearized operator around such element showing that, under a 
control on how far is the approximate state from being an exact stationary solution, 
a metastable behavior appears. 

This paper is organized as follows. 

In Section |2] we propose a construction for the family {?7^, V^} in the case of the 
Jin-Xin system. Then, we write the solution as 

u{x,t) = u\x,t) + U'{x;^{t)) 
v{x,t) = v^{x,t) + V'{x;^{t)) 

Hence we use as new coordinates the position of the shock layer ^, and the pertur- 
bation Y = {u^,v^). The couple {C,,Y) turns to solve an ODE-PDE coupled system 
of equations. Our approach consists in studying and approximation of such system, 
obtained by linearizing with respect to Y and by keeping the nonlinear dependence 
on ^, so that the o(y)-terms can be neglected. 

In Section [3| we analyze spectral proprieties of the linear operator £| arising from 
the linearization: we show that the spectrum of £^ can be decomposed into three 
parts: the first eigenvalue is real, negative and = 0{e~'"^^), C > 0, hence small 
as e — )■ 0"*'; all the others real eigenvalues are of order —C/e; all the remaining 
eigenvalues are complex with real and imaginary part less than —C/e, C > 0. Such 
estimates show that all of the components relative to all the eigenvectors except 
the first one have a very fast decay for small e, so that a slow motion occurs as a 
consequence of the size of the first eigenvalue. 

In Section |4| following the idea of [21] , we give an precise asymptotic expression 
for the first eigenvalue Xf^ of showing that it is exponentially small for e ^ 0. 

Finally, in Section [5} by using the spectral analysis performed in the previous 
Sections, we analyze the system for the couple {^,Y)] our main result is Theorem 
5.2, where we prove an estimate for the L^-norm of the perturbation Y, which can 



be used to decouple the (^, Y) system. 

Subsequently, we state a result on the behavior of the shock layer position, showing 
that it drifts towards the equilibrium location at a speed rate that becomes smaller 
as e — )■ 0. 

Finally, we numerically compute the position of the shock layer ^ at various time, 
showing that numerical results agree with the analytical results. 
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2. General Framework 



Let us consider the Jin-Xin system 
' dtu + dxV = 

1 



(2.1) 



dtv + a^d^u = ^(/(w) - V 
u{±i, t) = u± 

u{x, 0) = Uq{x 



v{x,0) = vo{x) = f{uo{x)) 



xe I, t > 

t > 

X e I 



for some flux function / chosen so that assumptions ( |1.3 ) hold. System (2.1) can 
be rewritten as 



(2.2) 
where 



dtZ = 7'\Z\ Z 



t=o 



r\z\ :- 



Vl[Z] 

n[z] 



- d,r.v 



-o?d^u + -{f{u) - v) 



We are interested in studying the behavior of the solution to ( 2.2[ ) in the relaxation 
limit, i.e. e — )■ 0. We assume that there exists a one-parameter family of functions 
{f/"(jO,'t^'(-;0}ee/ such that {U'{-]C),V'{-]C)) = {U',V') for some ^ e /. When 
^ 7^ ^, an element of this family can be seen as an approximate stationary solution 
to the problem, i.e. J-'^[f/^(-; 0^ ^"^(s 0] — as e — )■ in an appropriate sense to be 
specified. Moreover, we require that, in the relaxation limit, {U^{-] ^),V'^{-] ^)) — )■ 
(^hyp(-; 0' ^yp('5 0)- We remark that, once the one-parameter family of functions 
{U^i'^O} is chosen, the couple {U^{-; C,),V'^{-] C,)} is univocally determined by the 
relation 

V = -ea^d^U' + f{U') 



Example 2.1. In the case of Burgers flux, i.e. f{u 
(2.1) satisfies 

(2.3) ea^dxU = — ^, v 



a stationary solution to 



=PM*, for some u* > 0. An approximate solution 



with boundary conditions u{±l 
f/^(a;;^) to the first equation of (2.3) is obtained by matching two different steady 
states satisfying, respectively, the left and the right boundary conditions together 

^ = (see [201 Example 2.1]). In formula 

fc_tanh(A;_(^-a;)/2£) in(-£,0 

k+ tanh {k+{^ - x) /2e) in (^, i) 

1, and k± are chosen so that the boundary conditions are satisfied 



with the request f/*^ 
(2.4) U'ix;0 



where a 
(2.5) 



A;±tanh(^^(eT 
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Moreover, by the condition v = Hy, we have 



e/2 in(-£,0 
kl/2 in(e,£) 



2.1. The hnearized problem. As aheady stated before, in order to describe the 
dynamics generated by an initial configuration localized far from the steady state 
{U^, V'^), we assume to have a one-parameter family 

W%x;m) ■= {U'ix;m),V%x;mhei 

parametrized by ^ ^ I, such that the couple (f/^(x; ^(t)), ^(t))) is an ap- 
proximate stationary solution to (2.1), in the sense that it satisfies the stationary 
equation up to an error that is small in e. More precisely, following the idea firstly 
introduced in [2^, we assume that there exist two families of smooth positive func- 
tions fll = flK^) and f2| = uniformly convergent to zero as e — )■ 0, such that, 
for any ^ € /, the following estimates hold 



(2.6) 



m-),vnw^i;0])\<mM 
m-),n[w^i;0])\<mM\ 



VV^ G C(/) 

y^j G c{i) 



Once a one-parameter family {W^(-;,^)} satisfying (2.6) is chosen, we look for a 
solution to (2.1) in the form 

Thus we describe the dynamics in a neighborhood of the family {U'^{-; C,),V'^{-] ^)} 
using as coordinates the parameter ^ and a distance vector Y = {u^, v^) , determined 
by the difference between the solution (u, v) and an element of the approximate 
family. Substituting in (2.1), we obtain 

dtu' + d^U'i-; 0^ + d.V'i-, + dy = 
Since f{U' + u^) = f{U') + f'{U')u^ + o{u^) , we get 



(2.7) 
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where 



Vt[W] := -d^V 



Q'[u] ■= o{u) 



Example 2.2. Let us recall the Example 2A_ where we construct an approximate 
stationary solution for the Jin-Xin system with f{u) = v? /2 and a = 1. We compute 
in this specific case 

From the explicit formula for U^{x;^) given in ( 2.4[ ) we get, 'P|[W'^] = 0. 

On the other hand, —dxV^{x;^) = ed^xU^^ix;^) — dxf{U'^{x;^)). By direct sub- 
stitution, we obtain the identity 

in the sense of distributions. We also have 

In order to determine the behavior of P[[W^(-; ^)] for small e, we need an asymptotic 
description of the values k±. Following the idea of [20j, let us set k± := =F^±(1 + h±) 
and A± := £=f^. Relation (2.5) becomes 

tanh + 



2e ' ') l + /id 
Therefore, the values h± are both positive and then 

tanh f T-T^^ I < 



V 1e J ~ l + h± 
that gives the asymptotic representation 

(2.8) /i± < , , V , , - 1 = = 2e^"±^±/" + l.o.t 

^ ' tanh(^M±A±/2e) eT«±^±/^ - 1 

where l.o.t. denotes lower order terms. Finally 

1 

where u± = =Fu* for some u* > 0, so that we end up with 

(2.9) ld.U%=^ < ^(e-"*('^+«)/^ - e-"*('-«)/^) + l.o.t. 

showing that this term is exponentially small for e — )■ and it is null when ^ = 0, 
that corresponds to the equilibrium location of the shock when f{u) = u'^/2. 
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In this case, if we neglect the lower order terms, we can write an asymptotic 
formula for fif , that is 

(2.10) niiO ~ ^(e-"*(^+«)/" - e-"*(^-«)/") 



It follows that the hypothesis (2.6) is satisfied in the special case of f{u) = u /2. 

We can also numerically compute the limit of the solution (f/^, for e — t- 0+. 
For fixed ^, we observe that, as e becomes smaller, the transition between U- and 
becomes more sharp, while v tends to f{u*)6x=^, according to the fact that, in 
the limit e — )■ 0+, the solution (?7^(-; i), V^{-; ^)) converges to {Uhyp{-;^), Hiyp("; 0)- 
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Figure 2. Profile of the stationary solution {u,v) = {U^,V^) when 
f{u) = v? /2. The steepening of the shock layer and the convergence to 
a Delta function of u as e becames smaller are depicted. 



12 



M. STRANI 



Let US go back to the system (2.7). From now on, we set {u^,v^) = {u,v) and 



Y 



-a^d,u+-{f{U')u-v) 



Moreover, we introduce the following notation: if ip, (p E C, then {ip,4>) := Jj'ip4>, 
while if if) = (V^i,V'2) and (f) = (0i,02), then (t/?, 0) := + {ip2,h)- 

Miming the approach of pOj let us assume that, for any ^, the linear operator £| 
has a sequence of eigenvalues A| = A|(^) with corresponding (right) eigenfunctions 
<Pl = <pl{^,-) (for more details see Section 3). Denoting by i/^l = i/^KC,,-) the 
eigenfunctions of the corresponding adjoint operator and setting = Yk{C,]t) := 
(i/'K-; ^), y(-, t)), we impose that the component Yi = {ui,vi) is identically zero. 
More precisely, since we will prove that the first eigenvalue is real, negative and tends 
to zero as e — )■ , we need to solve the equation in a subspace in which the operator 
doesn't vanish. To this aim, we set an algebraic condition ensuring orthogonality 
between tj^l and Y, in order to remove the singular part of the operator Thus, 
denoting by Iq the initial datum of the perturbation, we have 

(2.11) ^(^l(..^(t)),Y{;t)) = and (V.^,(.; ^o), Fo(-)) = 



dt 



so that 



m;0,dtY) + {d^il)li;0%Y) 







Since i/jI = (ip^jip^) is the first (left) eigenfunction, there holds = ■^I'^'i) that 

is 



V 



Ai 



Hence, from (2.11) we get 
/ 



+ 



.of + n[w^(-,0] 



Since {il)\,C\Y) = XfiiklY) = 0, we have 

and we end up with a scalar differential equation for the variable ^, that is 

_ W(-;0,Q1^]+^l[w-(-,0]) + (V;r(-,0,Pi1w-(-;0]) 

dt a'^{C,,u,v) 



(2.12) 
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where 



Since we are interested in the regime y ~ 0, the equation (2.12) is approximately 
solved for small Y. Thus the term u,v) is expanded for u,f ~ 0, yielding 



1 



1 {d^i^l{-,0,Y) 



where 



Now, for sake of simphcity, let us call o.q{0 := {il^\{-;i),d^yV^{-,S,)). Thus we end 
up with the nonlinear equation for ^(t), which reads 



(2.13) 
where 



(2.14) 



p 

Ri : 



{d^rA-;0,Y)' 



[c^m-{dm-^o,Y)]{amy 

i^lQ^Y]) 



[Q'[Y] = {0, Q^u]) 



Equation (2.13) has to be coupled with the equation for the perturbation Y. To this 
end, (2.7) is rewritten in the form 

(2.15) 



dtY = qY - d^W{-, 0^ + -^'[W^] + Q'[Y] 



Using (2.13), we end up with the following equation 

(2.16) dtY = {q + Ml)Y + H'{x;0+T^'[Y,^] 
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MIY 



-d^V%--Op'[tY] + Q^[u] 

Hence wc obtain the following coupled system for the shock layer location ^{t) and 
the perturbation Y 



(2.17) 



dt 



r(0 1 + 



dtY = {q + Ml)Y + H%x; + n'[i, Y] 
Example 2.3. Let us consider the Jin-Xin system, for which one obtains 

For what concerns the linear operator, setting b^[x]^) := f'{U^{-,^)), we get the 
following expressions 



qY :-- 



e 



Cl'*Y :-- 



\ 



e 

dxU V 

e 



complemented with Dirichlet boundary conditions. To obtain an asymptotic expres- 
sion for the function 6'^(^), we need to approximately compute the functions i/'i and 
9^W^. As usual, we refer to the case f{u) = /2. 

For £ ~ 0, the function is close to the eigenfunction 'j/j" = (V'l'") V'l'") of the 
operator relative to the eigenvalue A = 0, with 

For example, in (— £, ^) we have 



( „2 



i! o ,0,1; I ^— jO.v r\ 

e 



e 

L = 0, = 



1 



■ 0,-!; 







That is = A{1 



(x+xo)\ 



) and ip'l'^ = edxi^i'^- By imposing the conditions on 



;0,U 



the boundary and on the jump, and by doing the same computations in the interval 
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i) we obtain 

(1 - e«+(^-«)/«'^)(l - e-«-(^+^)/«'^) X < e 
(1 - e-"-(^+«)/"'^)(l - e«+(^-^)/«'^) x>^ 



so that = ('01 , V'l) ~ (1)0) for e ~ 0. Furthermore, with the approximation 
U%x;^) ~ Uhyp{x;^) and ~ Vhyp{x), we have 

/^MX(,.e+.)(a;) 

V%x;^ + h)-V%x;0 1 


so that S^?/"^ and d^V^ converge to — [wj^g and — respectively as e — )■ in 
the sense of distributions. Thus, since {i/)l,d^'W^) ~ — {u}, we deduce an asymptotic 
expression for the function 6'^ 



With the choice of = (f/^, V^) proposed in Example 2.1 , such expression becomes 
(2.18) r(0 ~ ^(e-"*('+«)/^ - e-"*('-«)/^) 



3. Spectral analysis 

In this Section we analyze the spectrum of the linearized operator £| in order to 
determine a precise description of the location of the eigenvalues. 
We recall the expression of the operator 

(— dxV 

so that the eigenvalue problem = A$ reads 

1 

X^ = -a'd,<f + -{nun^-^) 

complemented with Dirichlet boundary conditions. Hence, by differentiating the 
second equation with respect to x, we obtain 

(3.2) ea'^d^^if - d^fiW)^) = A(l + eX)^ 
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Then we are interested in studying the eigenvalue problem for the differential linear 
diffusion-transport operator 



(3.3) 

In [20] it is proven that, under opportune hypotheses on the behavior of the function 
b'^{x; ^) in the limit e — )■ 0, the eigenvalues of C'^^'"^'^ have the following distribution 

-Ce-'^/" < < and XT" < -j VA; > 2 

More precisely, the following Propositions are proven in j20j (for more details, see 
Proposition 4.1 and Proposition 4.3). 

Proposition 3.1. Let he a family of functions satisfying the assumption: 
AO. There exists C > 0, independent on £ > 0, such that 



dx 



< C 



If there exist ^ G 6+ < < 6_ and C > for which \¥ — < Ce, where 

Ip is the step function jumping from to b'^ , then there exist constants C,c > 
such that -Ce-"l^ < Af ^ < 0. 

Proposition 3.2. Let b^ he a family of functions satisfying the assumptions: 
Al. ¥ E C°[— £, £], b^ is twice differentiahle at any x ^ ^ and 

d%' , „ , dV d%' 

— ,— <0<¥ m(-£,0, and 6^ — < < — in(e,£) 

A2. For any C > there exists cq > such that, for any x satisfying |x — ^| > 
cqS, there holds 

dV 

|6^-6°|<C£ and e 



dx 



< C 



A3. There exist the left /right first order derivatives ofV at ^ and 

\dV 



lim inf e 

e-s>0+ 



dx 



> 



Then there exist a constant C > such that X^" < —C/e for all e sufficiently small. 
Remark 3.3. When f{u) = m^/2, then ¥{x;^) = f/^(x;^). With the choice of [/^ 



proposed in Example 2.1 , we can easily check that hypotheses AO-1-2-3 are verified. 



From (3.2), we observe that A is an eigenvalue of if and only if A'"*^ := A(l + eA) 
is an eigenvalue for the operator C'^''"^'^ defined in (3.3). Hence, if A 



Xl^ is an 



eigenvalue of then there exists an eigenvalue XT'" such that 
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SO that 



(3.4) Ag = -^±^Vl + 4eAr 

Hence, if AJ^*"^ > then X:^^ G M. Moreover, since A^***^ are negative for all n G N 

9 \vsc 9 xi^sc 

(3.5) A^;:^ = < 0, A;:^ = , < 



Thanks to Propositions 3.1 and 3.2| we know that A^'^'^ > S'lid A^**^ ~ e '"/^ as 



e — 7- 0. Thus, from (3.4) and (3.5), there exists a constant C such that 

-e-^'/e < < A(^ < 

Moreover, if for some n > 1 there exist others A^'*'^ such that A^*'^ > then they 
are of order 1/e, so that 

AS < -C"/e 

On the other hand, if X^^'^ < then A;(^ G C. More precisely 

Proposition 3.2 assures that there exists k > 2 such that X'"^'"^ < n > k, 

so that -Re(AS) and /m(A;:[^) are terms of order 1/e. For example, if = 2 and 
we take into account Aj'"^ < 0, the corresponding eigenvalues for £| verifies 



^e(A2^,l) = --, Im{XiX) = ±^V|1 + 4eA^ 
Moreover, for A3^, since |A3''^| > |A2*^|, we have 

Re{Xi^^) = ~= Re{Xil) \Im{\i%)\ = > |/m(Ag)| 

Figure |3] shows the connection between the two spectra when A; = 2, so that only 
the first two eigenvalues of £^ belong to M. 
Hence, the following proposition holds 

Proposition 3.4. Let b"^ be a family of functions satisfying assumptions AO-1-2-3 
for some ^ G (£, t) and for some 6+ < < 6_ . Then the spectrum of the linearized 
operator £| can be decomposed as follow 

1. Xi^ G M and -e"^'/" < A(^ < 

2. Xi^_ G M and Xf^ < -1/e 

3. There exists k > such that 

Xi^^eR and Xi%<-C"/e, V n = 2, 1 + A; 
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C 




1 

\jx 












At 

■^t 1 

1 












1 

1, JX I ' 

-^3 ■•■ 

1 











Figure 3. The spectra of the operators and L^''"^'^. 
4- Xnj, eC for alln>2 + k and 

Remark 3.5. In [12], Kreiss G. and Kreiss H. performed the spectral analysis for 
the operator 

CeU := ed^:y,u - d^{f'{U'{x))u) 
arising from the linearization around the exact steady state lJ^{x) of 

dtu = ed^^u - d^J{u) 

proving that all the eigenvalues are real and negative. By using this result and our 
spectral analysis, if we linearize the system ( |2.1 ) around the exact stationary solution 



{U^^V^), we can prove that the real part of all the eigenvalues of the linearized 
operator is negative, so that the steady state (t/^, V^) is asymptotically stable with 
exponential rate. 

4. Asymptotic estimates for the first eigenvalue 
In this Section we want to study the behavior in e of the principal eigenvalue of the 



operator £^ associated to the linearization of (2.1 ) around an approximate stationary 
solution. Since usually the metastable behavior is the result of the presence of a 
first small eigenvalue, our aim is to determine an asymptotic expression for Af^. We 
have already emphasized the fact that A"^^ is an eigenvalue of the nonlinear Jin-Xin 
system if and only if A"**^ = A"^'^(1 + £:A'^^) is an eigenvalue for the following operator 

(4.1) C''"'' := ed^:,u - d^{f{U%x-i))u) 
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where U'^{x;C,{t)) is an approximate stationary solution for the scalar conservation 
law 

dtu = ed^^u - d^f{u) 
u{±i,t) = ^u*, u{x,0)=uo{x) 
In particular 



(4.2) 



(4.3) 



JX\ 



2|Af 



1 + JI + AeXT 



In [20], in the special case f{u) = f/^(a;; .^(t)) is defined in (2.4). In that paper 
it is proven that, for e ~ 

* 2 



2e 



so that 



(4.4) 



AS (6 1 



+ e 



1 + VI - 2m*2 [e 



This formula shows that the principal eigenvalue of the Jin-Xin system with f{u) = 
M^/2 is exponentially small in e. 

In order to determine an asymptotic expression of the first eigenvalue of the 



operator (4.1) for a general class of fiux function f{u), we refer to the paper of 
Reyna L.G. and Ward M.J., [23]; here the authors use the method of matched 
asymptotic expansions (MMAE) to determine an approximate stationary solution 
to dO). 



Miming their approach and computing the same calculations of [24j, with the 
appropriate changes due to the fact that the study of our equation is performed in 
the interval (— £, i) instead of (0, 1), the leading order MMAE solution for e — )■ 0^ is 
given by a function Us{x; ^) ~ Us[e~^{x ~^)], where ^ G (— £, i) and the shock profile 
uJz) satisfies 



u',iz) = fMz))-fiu*), 



(4.5) 



uJz) u — a_e 



Usiz) ~ —u + a+e 



z 



-oo < 2; < oo 
-oo 



z 



-oo 



The positive constant v± and a± describe the tail behavior of Us{z) and are defined 
by 

^± = Tf'iTu*) 

± 



drj 



In particular, when f{u) = u /2, Us{z) = — m* tanh('u*z/2), according to (1.9). No- 



tice that the MMAE solution satisfies exactly the equation, while the boundary 
conditions are satisfy within exponentially small terms. Instead, the construction 
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presented in [201 Example 2.1] and followed in this paper gives a function [/^(x;^) 
that verifies exactly the boundary conditions and solves approximately the station- 
ary equation. 

The eigenvalue problem associated to the linearization around Ug is given by 

L0 = e^d^^cj) - V[e-\x - O]0 = A0 




(4.6) 



V{z) 



\[fMz))f + \rMz))u',{z) 



In pi] it is proven that the first eigenvalue of (4.6) has the following asymptotic 
representation (for details see [211 Formula (2.14)]) 



(4.7) 



Ar(o 



1 

2u* 



Finally, from (4.3), we get 
(4.8) |AS(OI 



+ a_z/_e" 



1 + ,/l - ^ [a+z/2e-.+.-K^-0 + a_z/2e---=" 



H^+«)] 



This formula shows that A/^ is exponentially small as e — )■ 0. We remark that, 
when f{u) = m^/2, = a_ = 2u* and z/_|_ = z/_ 



u*, so that (4.8) is the same as 



(4.4). 



5. The behavior of the shock layer position 



Let us consider the system (2.17) for the couple (^, Y) and let us neglect the o{Y) 
terms 



(5.1) 



^(0 1 + 



Yt = {q + MI)Y + H%x;0 



This system is obtained by linearizing with respect to Y and by keeping the nonlinear 
dependence on ^, in order to describe the slow motion of the shock layer position 
far from the equilibrium location ^. 

We complement the so called quasi-linearized system (5.1) with initial data 

(5.2) ^(0) = ^0 G (-^,^) and Y{x,0) = {uoix),vo{x)), uo,vo E L^i-i, i) 



The aim of this Section is to analyze the behavior of the solution to (5.1) in the 
limit of small e. Subsequently, we will prove a result that characterizes the behavior 
of the shock layer location, proving that it moves towards the unique stationary 
solution with exponentially small rate. 



(5.4) lim|l]^|^^ =0, lim|fi||^^=0 
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Before stating our result, let us recall the assumptions. Let the family {W^(-; ^)} 
be such that there exist two families of smooth positive functions fif and f2| such 
that 

m-),n['^%;0])\ < mm,^ g c{i) 
^ ' ' m-),n[w^i;0])\ < mm.^ vv^ g 

We also assume that is asymptotically a solution, i.e. we require that 

(5.4) lim 

uniformly with respect to ^. 

Example 2.2| show that (5.3) and (5.4) are verified in the case of the quadratic 
flux f{u) =lF/2. 

For what concern the eigenvalues of the linear operator we have proven that 
there exist two positive constants Ci,C2 independent on ^ such that 

(5.5) Ag(0-i?e[Ag(0]>Ci, -e"^^/^ < Ag(0 < V^G (-£,£) 
Additionally, we assume that there exists a constant C3 > such that 

(5.6) QKe) + fil(0<C3|Ag(OI, ^ie{-lA) 

By comparing the asymptotic expression for Af^ given in (4.4) with the one for VL\ 
and Vl\ obtained in Example 2.3, we can easily check that hypothesis (5.6) is verified 
for the Jin-Xin system when f{u) = /2. 

Finally, concerning the solution Z = {z,w)'^ to the linear problem dtZ = C^^Z, we 
require that there exists z/ > such that for all ^ G (—£,£), there exist constants 
and C such that 

(5.7) |(z,^)(t)L <Q|(zo,^o)L.e-^*, Q < C V^G (-£,£) 



Remark 5.1. The assumption that < C for all ^ means that the estimate (5.7) 
holds uniformly in ^. Since ^ belongs to a bounded interval of the real line, if we 
suppose that ^ (-)■ is a continuous function, then there exists a maximum C in 

5.1. Estimate for the perturbation Y . Our first aim is to obtain an estimate 
the perturbation Y . We recall that 

(5.8) dtY = {C\ + M\)Y + H'{x-i) 

where 



^i1w^(-;0]-55f/^(-;er(0 



MiY 



In particular, is a bounded operator, such that 

(5.9) ||A^P£(L-M^) < C\9'm < Cmi^ + \ni\^^), Ve G i-i,^ 
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while 



)^l 



(5.10) \Hx^ <c^m^^+c2m 

For the special case of f{u) = both and are bounded by terms that 

are exponentially small in e, while, for a general class of flux functions f{u) that 
verify (1.3), the hypotheses we required assure that all the terms in the equations 



for the perturbation Y are small in e. 



Theorem 5.2. Let hypotheses (5.3)-(5.4)-(5.6)-(5.7) be satisfied. Then, fore suffi- 



ciently small, there exists a time T > such that, for all t < T, the solution Y to 



(5.8) is such that 
\Y\ 



for some positive constants Ci, C2 and 



.1(1 



0|L2 



snpXi'liO-Cmi 



4- 10^1 



> 



Proof. Since the operator £| + is a linear operator that depends on time, to 
obtain rigorous estimates on the solution Y, we need to use the theory of stable fam- 
ilies of generators, that is a generalization of the theory of semigroups for evolution 
systems of the form dtu = Lu. We will use some results of [22], which have been 
summarized in the Appendix A. More precisely, we want to show that £| + A^^ is 
the infinitesimal generator of a Co semigroup T^{t,s). 

To this aim, concerning the eigenvalues of the linear operator we know that 
■^^,+ (0 is negative and behaves like e~^/^ for all ^ G {—£, £), so that Af := sup^ ^i^iO 
is such that —e~^^^ < A^ < 0, and this estimate is independent on t. Hence, by 
using Definition 6.1 and Remark 6.2 (see Appendix A), we know that, for t G [0, T], 



since (|5.7l) holds, we get 



is the infinitesimal generator of a Cq semigroup iS^(f)(s), s > 0. Furthermore, 



so that the family 

u = — |A^|. Furthermore, since 



i^£) is stable with stability constants M = C and 



Theorem 



6.3 



\\M 

(see Appendix A) states that the family {£ 



stable with M = C and a; = -|Af | + C{\ni\ 



m 
< 0. 



From Theorem 6.8 (see Appendix A), we can define 7^(t, s) as the evolution system 
of dtY = {CI + Ml)Y, so that 

(5.11) Y{t) = ri^{t,s)Yo+ I r^{t,r)H%x;^{r))dr, 0<s<t<T 

J s 

Moreover, there holds 

mt,s)\\ < Ce'^'^'-'\ := \^\\~c{\n\l^ + 



> 
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Finally, from the representation formula (5.11) with s = 0, it follows 



(5.12) 



dr, 0<t<T 







so that, by using (|5.10|), we end up with 
(5.13) 



Y\^,it) < [Ci\ni\^^+C2\ni\^^] (l-e-^^*) +e-^^*|Fo|L2 



□ 



Remark 5.3. In the special case of Burgers flux, /x^ is going to zero exponentially 
as e — 0, since behaves like e~^/^ and from the explicit formula of and fig in 
Example 2.2, In the general case, assumption (5.4) assures that /i^ — )■ as e — )■ 0. 



Remark 5.4. To apply Theorem 6.8, we need to check that the domain of £| + A^^ 
does not depend on time, and this is true since £^+A^| depends on time through the 
function U^{x; ^(t)), that does not appear in the higher order terms of the operator. 
More precisely, the principal part of the operator does not depend on ^(t). 

5.2. Slow motion of the shock layer. An immediate consequence of the estimate 
(5.13) is that, for \Y\ ^ < M for some M > 0, the function ^(t) satisfies 



^ = r(0(l + r) with \r\<Ci\Ql\^^+C2\Q: 
More precisely, we can prove the following 



2\l° 



+ e 



Y. 



l2 



Proposition 5.5. Let hypotheses (5.3), (5.4), (5.6) and (5.7) be satisfied. Assume 
also 



(5.14) 



s e%s) < for any s G /, s ^ and ^'(0 < 



Then, for small e and |1o|^2j solution ^ converges to ^ as t +oo. 



Proof. Thanks to assumption (5.4), for e and |^o|^2 sufficiently small, estimate (5.13) 
holds globally in time. Hence, for any initial datum the location of the shock 
layer satisfies 

dz 



(5.15) 
where 



eHz) 



[1 + r{s))ds 



\r{t)\ < Ci\ni\^^+C2\ni\^^+e-'''\Yo\L2 



More precisely, if e and |^o|^2 sufficiently small, the dynamics of the shock 
location ^ is well described by 



(5.16) 



d^ 
dt 
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Therefore ^ converges to as t — )■ +00, and the convergence is exponential for any t 
under consideration, since, by means of the standard separation of variable method, 
we obtain the following estimate for the shock layer location 

(5.17) l^^t) - ^\ < \^,\e^'\ 

where 6'^'(^) -)■ as e -> 0. 

□ 



Formula (5.17) shows the slow motion of the shock layer for small e. Precisely, the 
evolution of the collocation of the shock towards the equilibrium position is much 
slower as e becomes smaller. 

For example, when f{u) = ^ = and ^^^'(0) ~ e~^/^ (see formula (2.18)). 

We also emphasize that hypothesis (5.14) are verified in the case of the Jin-Xin 
system with f{u) = 

The following table shows a numerical computation for the location of the shock 
layer for different values of the parameter e and f{u) = u'^/2. The initial datum for 
the function u is Uo{x) = |a;^ — x — |. We can see that the convergence to ,^ = is 
slower as e becomes smaller. 



The numerical location of the shock layer £^{t) for different values of the parameter e 



TIME t 


at), 8 = 0.1 


^t), e = 0.07 


£ = 0.055 


e = 0.04 


e = 0.02 


0.2 


-0.4008 


-0.4020 


-0.4029 


-0.4040 


-0.4059 


1 


-0.3314 


-0.3345 


-0.3360 


-0.3374 


-0.3389 


10 


-0.3070 


-0.3263 


-0.3304 


-0.3320 


-0.3326 


10^ 


-0.0103 


-0.1600 


-0.2562 


-0.3181 


-0.3325 


lO'^ 


-1.9725 * 10-1^ 


-0.0084 


-0.1115 


-0.2531 


-0.3320 


0.5* 10^ 


-1.9725 * 10-1^ 


-2.2102* 10-" 


-1.5057* 10-"* 


-0.0379 


-0.3099 



Figure |4] shows the dynamics of the shock layer (i.e the dynamics of the solution 
u to ( 2.1[ )), obtained numerically. When e = 0.1, the shock layer location converges 
to zero very fast: as we can also see from the table, when t = 10^, the value of ^{t) 
is already very close to zero. On the other hand, when e becomes smaller the shock 
layer location moves slower and it approaches the equilibrium location only for very 
large t. Finally, Figure [5] shows the profile of the shock layer for the fiux function 
f{u) = u'^/4:, that still verifies hypotheses (1.3). 



6. Appendix A 

In this section, we briefiy review some results on the theory of evolution systems 
by A. Pazy [221 Chapter 5]. For more details and for the proofs of the Theorems, 
see [22], Theorem 2.3, Theorem 3.1, Theorem 4.2. 

Let X be a Banach space. For every < t < T, let A{t) : D{A{t)) C X X he 
a linear operator in X and let f{t) be an X valued function. Let us consider the 
initial value problem 

(6.1) dtu = A{t)u + f{t), u{s)=uo 0<s<t<T 
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Figure 4. The shock layer profiles for different times and different values 
of the parameter e. 



In the special case where A{t) = A is independent of t, the solution to (6.1) can be 
represented via the formula of variations of constants 



u{t) = T(t)uQ + [ T{t- s)f{s) ds 
Jo 



where T{t) is the Cq semigroup generated by A. In [22] it is shown that a similar 
representation formula is true also when A{t) depends on time. 

Definition 6.1. Let X a Banach space. A family {A{t)}te[o,T] of infinitesimal gen- 
erators of Co semigroups on X is called stable if there are constants M > 1 and u 



26 



M. STRANI 



u 




.1 2 I ^ I ^ I ^ I ^ I I I 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

E=0.07, f(u)=u'*/4 

Figure 5. Profiles of the shock layer at different times with a convex flux 
function f{u). 

(called the stability constants) such that 

{cu, +00) C p{A{t)), for t e [0, T] 

and 

||nJ^ii?(A : A{tj))\\ < M{X-u)~\ 

for X > uj and for every finite sequence <ti < t2, tk < T , k — 1, 2, 

Remcirk 6.2. If for t e [0, T], A{t) is the infinitesimal generator of a Co semigroup 
St{s), s > satisfying ||5't(s)|| < e^^, then the family {A{t)}te[o,T] is clearly stable 
with constants M — 1 and cu. 
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The previous remark means that, if for every fixed t G [0,T] the operator A{t) 
generates a Cq semigroup St{s), and we can find an estimate for [[^((s)!! that is in- 
dependent of t, then the whole family {A{t)}t^[Q T] is stable in the sense of Definition 

Theorem 6.3. Let {A{t)}t^[Q^T] be a stable family of infinitesimal generators with 
stability constants M and u. Let B{t) , < t < T be a bounded linear operators 
on X . If < K for all t < T, then {A(t) + i?(t)}(g[o,r] is a stable family of 

infinitesimal generators with stability constants M and u + MK. 

In order to prove the existence of the the so called evolution system U{t, s) for 
the initial value problem (6.1 ), let us introduce X and Y Banach spaces with norms 
II II X, II II y respectively. Moreover, let us assume that F is a dense subspace of X 
and that there exists a constant C such that HtfUx < C'||w||y for all w 

Definition 6.4. Let A be the infinitesimal generator of a Co semigroup S{s), s > 0, 
on X . Y is called A-admissible if it is an invariant subspace of S{s), and the re- 
striction S{s) of S{s) to Y is a Co semigroup on Y . Moreover, A is the infinitesimal 
generator of the semigroup S{s) on Y , and it is called the part of A in Y . 

Next, let us fix t G [0,T], and let A{t) be the infinitesimal generator of a Cq 
semigroup St{s) on X. The following assumptions are made 

(HI) {A{t)}ti^[Q,T] is a stable family with stability constants M and u. 

(H2) Y is y4()f:)-admissible for t G [0,T] and the family {yl(t)}tg[o,T] is a stable 
family in Y with stability constants M, a). 

(H3) For t G [0, T], y C D{A{t)), A{t) is a bounded operator from Y into X and 
t —7- A{t) in continuous in the B{X, Y) norm. 

Remark 6.5. The assumption that the family {A{t)}ti^[Q,T] satisfies (H2) is not al- 
ways easy to check. A sufficient condition for (H2) which can be effectively checked 
in many applications states that (H2) holds if there is a family {Q{t)} of isomor- 
phism of Y onto X such that ||Q(t)||y^x and ||Q(t)^^||y^x are uniformly bounded 
and t ^ Q{t) is of bounded variation in the B{Y, X) norm. 

Remark 6.6. Condition (H3) can be replaced by the weaker condition 

(H3)' For t G [0,T], Y C D{A{t)) and A{t) G L\[Q,T]- B{Y, X)). 

Theorem 6.7. Let A{t), < t < T be the infinitesimal generator of a Cq semigroup 
Stis), s > on X . If the family {A(t)}(g[o,T] satisfies the conditions (H1)-(H3), 
then there exists a unique evolution system U{t, s), < s < t < T , in X satisfying 

(6.2) \\U{t, s)\\ < Me'^^^-'\ for < s < t < T 



Moreover, if f G C{[s,T], X), the solution to (6.1) can be written as 
(6.3) u{t) = U{t,s)uo+ [ U{t,r)f{r) dr 



28 



M. STRANI 



for alio <s<t< T. 

One special case in which the conditions of Theorem |6.7| can be easily checked is 
the case where the domain of the operator D{A{t)) = D is independent on t. In this 
case we can take D as the Banach space which we denote by Y, and the following 
Theorem holds 

Theorem 6.8. Let {A{t)}t^[Q^T] be a stable family of infinitesimal generators of Cq 
semigroups on X. If D[A{t)) = D is independent on t and for uq G D, A{t)uo is 
continuously differentiable in X, then there exists a unique evolution system U(t,s), 



< s < t < T , satisfying (6.2). Morevoer, if f G C{[s,T], X), then, for every 



Uq G X, the initial value problem (6.1) has a unique solution given by (6.3) 
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